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Abstract. We present the results of a comprehensive Spitzer survey of 70 

radio galaxies across 1 < z < 5.2. Using IRAC, IRS and MIPS imaging we de- 
termine the rest- frame AGN contribution to the stellar emission peak at 1.6/xm. 
The stellar luminosities are found to be consistent with that of a giant elliptical 
^ ■ with a stellar mass of 10^^~^^ Mq. The mean stellar mass remains constant at 

\jQ ' ^ lO^^*^ Mq up to z = 3 indicating that the upper end of the mass function 

fSJ , is already in place at redshift 3. The mid-IR luminosities imply bolometric IR 

CSJ ■ luminosities that would classify most sources as ULIRGs. The mid-IR to radio 

^^ ' luminosity generally correlate implying a common origin for these emissions. 

^D , The ratio is higher than that found for lower redshift, i.e. z < 1, radio galaxies. 

\o 

O 
I ' 1. Introduction 

o 

^ ' Classical radio galaxies are the quintessential type II active galactic nuclei (AGN) : 

C^ . accreting super-massive black holes that have their continuum emission in the 

UV/optical/NIR absorbed by dust, thus primarily giving them the appearance 
of normal star-forming galaxies at these wavelengths. The main evidence that 
^ I they host super-massive black holes comes from the high luminosities of their 

j^ . radio 'lobes' that are fed by radio jets originating from the host galactic nu- 

clei (Rees 1978). The lobe spatial extents (tens of kilo-parsecs) and the lu- 
minosities (ivi.4GHz ^ lO^^WHz"^) rule out emission from star-formation. To 
be more precise radio galaxies are Franhoff-Riley type II objects with edge- 
brightened radio lobes. In terms of the orientation unification scheme for AGN 
(Antonucci 1984) radio galaxies are analogous to radio loud quasars obscured in 
the UV/optical/NIR. 

Due to their large radio luminosities, radio galaxies were the predominant 
way to probe the distant universe until the advent of ultra-deep optical surveys 
in the last decade. In fact, radio galaxies were the first galaxies to be found 
above redshifts 1, 2, 3 and 4 (e.g., Stern & Spinrad). Since their first discovery 
it has been known that the optical hosts of luminous radio sources are primarily 
giant elliptical (gE and cD) galaxies (Matthews et al. 1964). In the more dis- 
tant universe, indirect evidence that this association remains intact comes from 
the detection of normal elliptical host galaxies with r^'^ law light profiles in 
-fTST/NICMOS observations of high-redshift radio galaxies (HzRGs) at l^z<2 
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(Pentericci et al. 2001; Zirni et al. 2003); the tendency for HzRGs to reside in 
moderately rich (proto-) cluster environments (Venemans et al. 2002; Stern et 
al. 2003); the spectacular (>100kpc) luminous Lya haloes seen around several 
sources, implying large gas reservoirs (Reuland et al. 2003; Villar-Martm et al. 
2003); sub-mm detections of HzRGs, implying violent star formation activity up 
to ~100MQyr-^ (Archibald et al. 2001; Reuland et al. 2004); and a few direct 
kinematic measurements of HzRGs (Dey & Spinrad 1996). The most compelling 
evidence of this association of HzRGs with the most massive systems, however, is 
the tight correlation of the observed near-infrared Hubble, oi K — z diagram for 
powerful radio sources (De Breuck et al. 2002, Rocca-Volmerange et al. 2004): 
HzRGs form a narrow redshift sequence which traces the envelope of radio-quiet 
galaxies and is well-modeled by the evolution of a stellar population formed at 
high redshift from a reservoir of 10^^ M©. With the more recent discovery that 
the stellar bulge and central black hole masses of galaxies are closely correlated, 
it is no longer a surprise that the parent galaxies of the most powerful radio 
sources occupy the upper end of the galaxy mass function (Maggorian 1998; 
Tremaine 2002). 

The peak of the stellar emission at 1.6 //m of elliptical galaxies has been 
found to be a reasonably robust measure of the stellar mass for old passively 
evolving systems. The mass-to-light ratio at this wavelength does not vary 
greatly for ages > IGyr. The Spitzer Space Telescope now allows us to observe 
this feature in distant sources in the rest frame. In particular the four bands 
of the IRAC instrument and the IRS 16 ^m peak-up imager straddle the rest 
frame 1.6 /im flux density at 1 <z< 5. 



1.1. Sample selection 

In order to investigate the formation and evolution of the most massive galaxies 
we have performed a Spitzer survey of 70 HzRGs in GO cycle 1. These HzRGs 
have also been carefully chosen to cover the full range of redshifts from z = 1 
to the redshift of the highest known radio galaxy (z = 5.2) and two orders of 
magnitude in radio luminosity, preferentially selecting targets with supporting 
data from HST, Chandra and SCUBA/MAMBO. By covering this parameter 
space, any trends with radio luminosity or redshift should be apparent. The 
observations consist of photometry in eight bands from the three instruments 
aboard Spitzer, exercising the full complement of imaging capabilities (IRAC: 
3.6, 4.5, 5.8 and 8/um; IRS: 16 /um; MIPS: 24, 70 and 160/um). Due to uncertainty 
in the ability of MIPS to image against the Galactic infrared background at the 
time of submission of GO Cycle 1, we chose only to image 26/70 HzRGs with 
MIPS. The other 54 sources have been applied to be observed in Spitzer Cycle 3. 
The IRS images are only for the 46 objects above z = 2 as below this redshift the 
8/im IRAC channel adequately covers the longward side of the 1.6 /xm bump. 

The 'SHizRaG' team keep track of this project through a private webpage 
which we intend to make public eventually. Currently a restricted version of the 
webpage is available here: 
|http : //spider . ipac . caltech . edu/staf f /seymour/SHizRaGs . html 
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Figure 1. (Left) Postage stamp images of 4C23.56 at z=2.48: the 4 IRAC 
bands (top row - wavelength increasing left to right) followed by the IRS 16 /im 
peak-up image and the MIPS bands (bottom row) showing a clear detection 
in each waveband out to 70 fim (Right) SED fitting of 4C23.56 Spitzer data 
using an elliptical galaxy template and black-bodies of various temperatures. 



Stellar Luminosities and Masses 



2.1. SED fitting 

In order to determine the 1.6 /im stellar luminosity, the contribution at this wave- 
length from hot, AGN heated dust needs to be ascertained. We have performed 
this analysis for the sample of 17 MIPS HzRGs which have 24 ^m detections as 
well IRAC observations. For the modeling presented here we have chosen to use 
just the Spitzer data. The principle reason is that shorter wavelengths, i.e. the 
rest-frame optical, are more likely to have significant contributions from young 
stellar populations, emission lines, and AGN continuum. 

We model the observed 3.6 to 24 /^m data using an elliptical stellar template 
and three black-bodies. The elliptical stellar SED is taken from the PEGASE 
(Fioc & Rocca-Volmerange 1997) software and corresponds to a passively evolv- 
ing 10^^ Mq stellar population formed at z = 10; i.e. for each redshift a template 
of the correct age was chosen although the SED does not evolve significantly af- 
ter 1 Gyr (i.e. z ~ 4.4). The three black-bodies where chosen to be at 60 K 
(analogous to the temperature of cold dust found from sub-mm observations), 
250 K and 600-1200 K. This last value was allowed to vary and the best tem- 
perature from the fitting applied. Full details of this modeling are presented in 
Seymour et al. {in prep.). 

Figure 1 illustrates the Spitzer observations and the SED modeling for one 
representative source, 4C23.56 at z = 2.48. The best fit hot dust temperature 
is 750 K and one can see that the 1.6 /im stellar peak has considerable AGN 
contamination. A more detailed analysis of the broad-band. X-ray to radio SED 
of 4C23.56 is presented in De Breuck et al. {in prep.). 

The other 53 HzRGs with no MIPS detections can at least provide upper 
limits to the stellar luminosities and, in some cases only the upper limits come 
from detections only in IRAC channels 1 and 2. So we fit a maximum elliptical 
template SED to the IRAC data. In some case the elliptical template fits quite 
well to IRAC channels 1-3 and in others the SED is steeply rising at the longer 
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Figure 2. (Top) HzRG Stellar luminosities plotted against redshift of each 
HzRG. The solid line indicates a of 10^^ Mq elliptical galaxy. (Bottom) 
Stellar masses from the SED fitting plotted against redshift. Most stellar 
luminosities indicate with stellar masses of 10^^ — 10^^ Mq. Stellar luminosities 
and masses derived from sources with MIPS 24/im detections are circled, the 
dots with out circles are 'nominal' masses from just IRAC data. Downward 
arrows indicate upper limits. 



wavelengths and the fit is restricted only by channel 1. In the former case the 
we calculate a 'nominal' stellar mass from the fit, but in the latter we may only 
derive upper limits. 



2.2. Results 

The derived rest-frame, AGN-subtracted, 1.6 ;um stellar luminosities are shown 
against redshift in Fig. 3 (top) for all our HzRGs. Also laid on the plot by 
a solid line is the stellar luminosity of a 10^^ Mq elliptical galaxy. The stellar 



luminosities imply stellar masses in the range 10 — 10 Mq with a mean mass 
of ~ lO"*^^'^ Mq (Fig. 3 bottom). This mean mass remains consistent out to z = 3 
(beyond which the parameter space becomes less well sampled) suggesting that 
the upper end of the mass function is already in place by at least z = 3. 
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Figure 3. Rest inid-IR luminosity against rest radio luminosity which ap- 
proximately correlate implying a common origin. This correlation may be 
explained by the AGN producing the radio jets and also heating the hot dust 
radiating in the mid-IR. The mid-IR/radio ratio increases at higher redshifts 
compared to low redshift samples (e.g., Ogle et al. 2006) most likely due to 
the hotter temperatures of the AGN- heated dust. 



3. Mid-IR luminosities 

Figure 3 shows the rest 6.75 /im luminosity against the rest 3 GHz luminosity 
for the 18 radio galaxies with MIPS observations. The wavelength of 6.75 /im 
was chosen as a fiducial mid-IR wavelength as it is the mean rest wavelength of 
the observed MIPS 24 /im band for our sample and also the wavelength of the 
LW2 filter from ISOCAM, allowing a direct comparison of derived luminosities. 
It is clear to first approximation that the two luminosities correlate and by 
implication have a common origin. This correlation makes sense if the radio 
luminosity comes from lobes induced by a jet from the AGN and the mid-IR 
comes from hot AGN- heated dust. 

The mid-IR luminosities are also all greater than W^^Lq implying bolomet- 
ric luminosities on ULIRG scales assuming local relations hold for these objects. 
These radio galaxies tend to have a higher mean mid-IR to radio luminosity ra- 
tio then those selected at lower redshift, e.g., Ogle et al. (2006) finds that z < 1 
radio galaxies tend to have -Zymid-mZ-Z^radio '^ 10 — 100. 



4. Conclusions and future ^vork 

We have presented a stellar-luminosity/redshift relation of HzRGs, a more phys- 
ical representation of the K — z diagram. This distribution seems to confirm the 
long held paradigm that radio galaxies are hosted by massive ellipticals out to 
high redshifts, and that the most massive galaxies are already in place by red- 



shift 4 and possibly earlier. We also observe a correlation between the infrared 
and radio luminosities which is unsurprising if they are both fueled by the AGN. 
Current on-going work includes mm/sub-mm observations with SCUBA, 
MAMBO and the CSO to constrain the cold dust component at longer wave- 
lengths and hence estimate the mass of this cold dust. Over-densities of sources 
around radio galaxies are being investigated to look for evidence of cluster for- 
mation (Zirm et al. in prep.). Over half the sources with 24 ^um images have 
over-densities of factors of 2-5 greater than that expected from 24 /xm source 
counts. These radio galaxies mainly lie at 1.5 < z < 2.5 where the strong 
6 — 8^m PAH feature passes though the 24 ^m MIPS band, enhancing the 
24;um flux density. 
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